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NOMENCLATURE 


c heat  content 

D hydraulic  conductivity  of  unsaturated  water 
Og  soil  constant  introduced  in  (16) 
h thickness  of  the  ice  iens 
L latent  heat 

p pressure 

q defined  by  (39) 

r defined  by  (39) 

5 entropy 

T temperature 

t time 

U„(|)  defined  in  (19) 

u s exp  ifiW),  defined  at  (1 7) 

V'  volume 

V flux  of  water 

W degree  of  saturation  by  movable  water 
w surcharge  on  the  ice  lens,  i.e.  frost-heaving  pressure 
X space  coordinate 

a thermal  diffusivity 

Ojo  thermal  diffusivity  of  water-saturated  soil 
fi  soil  constant  introduced  in  (16) 

0„(t)  defined  in  (42) 

K thermal  conductivity 

K20  thermal  conductivity  of  water-saturated  soil 

Xg  soil  constant  introduced  in  (40) 

p constant  introduced  in  (7)  and  formulated  in  (31) 

V soil  constant  introduced  in  (55) 

{ variable  defined  by  (20) 

p density 

X,  volume  occupied  by  both  soil  particles  and  unmovable  water  in  a unit  volume  of  the 
soil  mass 

Subscripts 
A air 

I In-situ  freezing 

i ice 

S segregation  freezing 

s soil 

water 

the  first  layer  in  Figure  3,  i.e.  the  ice  lens 
the  second  layer  in  Figure  3,  i.e.  the  unfrozen  soil 
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SEGREGATION-FREEZING  AS  THE  CAUSE  OF 
SUCTION  FORCE  FOR  ICE  LENS  FORMATION 


Shunsuke  Takagi 


INTRODUCTION 

The  most  enigmatic  problem  in  the  theory  of  frost 
heaving  is  the  generation  of  the  suction  force  that 
draws  pore  water  up  to  the  freezing  interface  of  a 
growing  ice  lens,  increasing  its  thickness  despite  the 
pressure  exerted  by  the  overlying  burden  of  the 
frozen  soil  and  surface  load  (see  cover  for  photograph 
of  ice  lens).  Most  of  the  current  literature  on  frost 
heaving  explains  the  suction  by  use  of  the  Laplace 
equation  in  capillary  theory,  which  gives  the  pressure 
difference  across  a curved  meniscus  boundary  of  two 
different  materials.  However,  capillary  theory  is  not 
yet  proven  to  be  valid  on  the  freezing  meniscus  of 
pore  water. 

In  the  case  of  a capillary  tube  containing  air  and 
water  separated  by  a meniscus,  the  molecules  com- 
posing the  meniscus  are  stationary,  and  the  Laplace 
equation  gives  the  pressure  difference  across  the 
meniscus.  The  molecules  composing  a freezing  menis- 
cus, however,  are  consuntly  renewed;  theory  does 
not  yet  prove  whether  the  Laplace  equation  is  valid 
or  not  on  such  a meniscus.  We  have  found  experi- 
mental evidence  indicating  that  the  Laplace  equation 
is  valid  on  a static  ice/water  meniscus  where  molecules 
are  stationary,  but  not  on  a freezing  meniscus  where 
molecules  are  renewed. 

In  their  definitive  theoretical  work  in  this  field, 
Everett  and  Haynes  (1965)  caution  that  their  theory 
of  ice  stress  derived  by  applying  capillary  theory  on 
ice/water  menisci  may  fail  when  kinematic  effects 
predominate.  Koopmans  and  Miller  (1966)  measured 
the  capillary  potential  of  the  ice/water  meniscus  and 
showed  that  the  resulting  curve  coincided  with  the 
soil  moisture  characteristics,  if  they  substituted  the 
ice/water  interfacial  tension  with  the  air/water  inter- 
facial tension.  They  took  24  hours  to  get  one  point 
of  daU.  Their  experiment  shows  that  capillary  theory 
applies  on  the  sutic  ice/water  interface.  Penner  (1967) 


and  Sutherland  and  Gaskin  (1973)  showed  that  the 
pressure  required  to  stop  ice  lens  growth  was  larger 
than  the  pressure  predicted  by  capillary  theory.  Their 
experiments  were  kinematic.  We  may  interpret  these 
studies  as  indicating  that  the  freezing  of  pore  water  is 
not  a sutic  effect  caused  by  the  capillary  pressure  but 
is  a kinematic  effect  caused  by  the  simulUneous  flows 
of  heat  and  water. 

We  have  developed  a concept  suitable  for  describing 
the  ice  lens  formation  by  using  the  theory  of  simul- 
Uncous  flows  of  heat  and  water  (Takagi  1959, 1963, 
1965, 1970, 1974, 1975, 1977),  which  is  systematized 
and  suted  in  this  report. 


SEGREGATION  FREEZING 

We  shall  introduce  segregation  freezing  as  the  agent 
for  creating  suction  force  to  draw  water  to  the  freezing 
front  and  exerting  frost-heaving  pressure  to  the  over- 
lying  burden. 

&)rte  (1962)  observed  that  ice  growing  upward  can 
carry  soil  particles  floating  on  the  surface  (Fig.  1 ).  The 
explanation  of  the  floating  of  a soil  particle  on  a 
heaving  ice  surface  is  possibie  only  by  assuming  that, 
between  the  surfaces  of  the  particle  and  the  ice,  there 
exists  a thin  layer  of  unfrozen  water  whose  molecules 
are  constantly  replenished  during  the  heaving  by  the 


Figure  1.  A particle  floating 
on  the  heaving  ice  surface. 


Flow  of 
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Figure  2.  Ice  lens  forming  on  the  thin  water  layer. 

influx  of  water  from  the  adjacent  reservoir  into  the 
freezing  front.  In  other  words,  we  should  recognize 
that,  adsorbed  or  absorbed  between  the  particle  and 
the  ice,  there  exists  a heterogeneous  layer  of  water 
whose  "thickness"  is  maintained  at  a certain  constant 
value  during  the  freezing  process.  It  should  be  stated, 
however,  that  theoretical  physics  cannot  yet  explain 
the  water  of  this  nature  and  we  are  quite  ignorant  of 
its  properties. 

In  Figures  1 and  2,  soil  particles  are  represented  by 
a rectangular  shape,  because  the  thickness  of  the 
heterogeneous  layer  is  clearly  shown  in  this  form.  In 
the  case  of  an  actual  more  complicated  shape,  the 
conceptual  correction  of  "thickness”  can  be  made 
easily. 

The  freezing  of  a thin  water  layer  generates  suction 
that  draws  water,  as  shown  in  Figure  1 , from  the  sur- 
rounding reservoir.  The  freezing  of  water  that  gen- 
erates suction  will  be  called  segregation  freezing.  In 
this  case  the  heterogeneous  water  adsorbed  or  absorbed 
between  the  particle  and  the  ice  freezes.  In  contrast, 
the  freezing  of  homogeneous  free  pore  water  will  be 
called  in-situ  freezing.  This  freezing  mechanism  does 
not  generate  suction;  i.e.  the  in-situ  freezing  front 
advances  with  the  progress  of  the  freezing. 

In  in-situ  freezing,  the  ice  pressure  and  the  water 
pressure  may  not  necessarily  be  equal,  but  mechanical 
equilibrium  is  established  between  ice  and  water.  Ice 
and  water  are  also  equithermal;  i.e.  they  are  in  thermal 
equilibrium.  Therefore,  the  equilibrium  of  in-situ 
freezing  is  twofold.  In  the  case  of  segregation  freezing, 
however,  thermal  equilibrium  is  established,  but  as 
shown  below,  mechanical  equilibrium  is  not.  Only  one 
type  of  equilibrium  is  present  in  this  case. 

The  in-situ  freezing  temperature  is  determined  by 
the  condition  that  the  three  phases  — ice,  water,  and 
vapor  — are  in  thermodynamic  equilibrium  (Takagi 
1959).  The  segregation  freezing  temperature  is  deter- 
mined, as  shown  in  the  following,  by  the  two-phase 
equilibrium. 


Let  us  consider  an  ice  lens  resting  on  soil  particles, 
as  shown  in  Figure  2.  If  the  uppermost  part  of  the  thin 
water  layer  freezes,  water  must  be  sucked  in  from  the 
neighboring  reservoir  to  recover  the  original  thickness 
of  the  thin  water  layer.  Then,  if  the  soil  particles  stay 
at  the  same  position  during  •'5e  freezing  process,  the 
surface  DD  rises  by  the  thickness  of  the  frozen  portion. 
This  is  our  explanation  of  frost  heaving.  According  to 
our  concept,  therefore,  an  ice  lens  grows  on  soil  parti- 
cles. 

Stress  in  the  thin  water  layer  that  sustains  the  weight 
of  the  ice  lens  plus  any  surcharge  on  it  acts  in  the  water 
layer  as  if  the  water  were  solid.  However,  to  calculate 
the  freezing  temperature  of  the  thin  water  layer,  i.e. 
the  segregation  freezing  temperature,  the  simpler 
thermodynamic  state  of  the  pore  water  underlying  the 
ice  meniscus  BMA  may  be  considered  instead  of  the 
complicated  thermodynamic  state  of  the  thin  water 
layer  under  the  flat  ice  surface  AB,  because  we  may 
assume  they  are  equithermal. 

The  pore  water  underlying  the  ice  meniscus  BMA 
and  the  ice  lens  overlying  BMA  are,  in  turn,  equithermal, 
but  not  in  mechanical  equilibrium.  The  weight  of  the 
ice  lens  is  not  supported  by  the  pore  water,  but  by  the 
thin  water  layers  and  the  soil  particles  underlying  the 
layers.  The  stress  of  the  ice  lens,  therefore,  is  inde- 
pendent of  the  pore  water  pressure.  In  soil  mechanics 
terminology,  the  pore  water  pressure  is  neutral  to 
mechanical  effects,  but  the  stresses  of  the  thin  water 
layers  and  the  soil  particles  are  effective.  They  belong 
to  different  categories  in  terms  of  mechanical  effects 
(Terzaghi  1942). 

For  simpler  treatment  we  replace  the  ice  stress  with 
the  ice  pressure  P\,  Then  we  can  describe  the  thermo- 
dynamic equilibrium  between  the  ice  lens  and  the  pore 
water  by  use  of  a formula  of  classical  thermodynamics: 

V^dP^-S^dT  = V^dP^-S^dT  (1) 

where  V is  the  specific  volume,  S the  specific  entropy, 

P the  pressure,  and  T the  temperature.  The  suffixes  w 
and  i refer  to  water  and  ice,  respectively.  Note  that  P^ 
and  P^  may  not  necessarily  be  equal  in  this  equation. 
(See  Takagi  (1965)  for  a treatment  dealing  with  the 
tensorial  ice  stress.) 

The  meaning  of  the  variations  denoted  by  the  total 
differentials  in  this  equation  must  be  clarified.  We 
choose  the  datum  state  (i.e.  the  starting  point  of  the 
variation)  to  be  the  state  of  in-situ  freezing.  The  tem- 
perature at  the  datum  state,  therefore,  is  the  in-situ 
freezing  temperature  Ti.  We  raise  the  pressure  of  the 
ice  at  the  datum  state  by 

dP\  = w + P|A  (2) 


2 
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where  w is  the  surcharge  overlying  the  ice  lens,  h the 
thickness  of  the  ice  lens,  and  pj  the  density  of  ice. 

We  do  not  change  the  pressure  of  the  pore  water: 

dP„  = 0.  (3) 

We  assume  that  the  soil  column  underlying  the  soil 
particles  is  incompressible,  so  that  no  disturbance  can 
intrude  into  the  system  during  the  proposed  pressure 
increase.  During  the  process,  wc  mainuin  the  thermo- 
dynamic equilibrium  between  the  pore  water  and  the 
ice  lens  by  keeping  (1)  valid,  but  leave  the  tempera- 
ture free  to  change.  Note  that  w + pj/t  is  the  surcharge 
on  the  ice  at  the  freezing  front,  which  may  be  inter- 
preted, if  frost  heave  actually  takes  place,  to  be  the 
frost-heaving  pressure. 

Thus  we  can  reach  the  final  stage  of  the  formula- 
tion. We  find  the  segregation-freezing  temperature 

Ti  = T^[^ -(w  + p.p)l(p^L)]  (4) 

by  letting  dr  = T^-Tf  andSj  = L/T^  in  (1),  where 
r,  is  the  in-situ  freezing  temperature  and  L the  latent 
heat.  Therefore,  T",  is  always  less  than  T^,  the  differ- 
ence being  determined  by  the  ice  pressure  increment, 
i.e.  the  frost-heaving  pressure. 

If  we  consider  that  the  stress  in  the  ice  lens  is  deter- 
mined by  the  configuration  of  the  ice  surface,  the 
stress  is  not  necessarily  uniform  in  the  ice  lens.  The 
nonuniform  stre«  caused  by  the  capillary  force  is  con- 
sidered by  Everett  and  Haynes  (1965).  The  difference 
of  the  capillary  forces  between  the  top  and  bottom 
menisci  is  considered  by  Loch  and  Miller  (1975)  to 
explain  the  cause  of  the  flow  of  ice  molecules  in  the 
growing  ice  lens.  However,  the  capillary  force  does 
not  seem  to  be  directly  related  to  the  crystal  growth. 

In  supercooled  water,  ice  crystals  grow  with  sharp 
edges  (Hobbs  1974,  Glen  1974)  and  frequently  form 
dendrites.  They  grow  against  the  chemical  potential 
gradient  in  the  solid;  their  growth  rate  is  determined 
by  the  heat  transfer  and  the  availability  of  the  growth 
material  in  the  liquid.  When  the  water  temperature 
is  very  close  to  the  ice  temperature,  however,  ice  grows 
into  the  water  forming  a smooth  ice  surface  (Glen 
1974).  The  growth  rate  in  this  case  is  still  determined, 
we  believe,  by  the  heat  transfer  and  the  availability  of 
the  growth  material  in  the  liquid,  although  we  could 
not  find  any  reference  that  clearly  sutes  this.  The  ice 
stress  caused  by  the  ice/water  menisci  does  not  seem 
to  be  a cause  of  crystal  growth. 

We  showed  (Takagi  1965)  that  the  ice  stress  given 
on  the  right-hand  side  of  (2)  is  the  normal  stress  com- 
ponent in  the  vertical  direction.  This  normal  stress 
may  be  interpreted  to  be  the  overall  representative 


value  of  the  ice  stress  in  the  segregation  freezing,  in  the 
same  sense  as  the  overall  representative  ice  stress  in  the 
in-situ  freezing  is  interpreted  (Takagi  1 959)  to  be 
atmospheric.  Obviously,  the  stress  of  the  ice  forming 
inside  a pore  of  soil  is  higher  than  atmospheric  pressure 
by  the  amount  of  the  capillary  pressure  caused  by  the 
curved  ice  surface.  The  formula  of  the  in-situ  freezing 
temperature,  derived  by  assuming  the  ice  pressure  to  be 
equal  to  atmospheric  pressure  is,  however,  confirmed 
experimentolly  (Schofield  1935,  Williams  1964,  Low  et 
al.  1968).  Atmospheric  ice  pressure,  therefore,  may  be 
the  overall  representative  value  of  the  internal  stress  of 
the  ice  freezing  in-situ,  and  choosing  atmospheric  pressure 
is  probably  a convenient  way  of  avoiding  the  variability 
of  the  internal  ice  pressure  in  in-situ  freezing.  Choosing 
the  ice  stress  expressed  by  the  right-hand  side  of  (2) 
in  the  formulation  of  the  segregation  freezing  tempera- 
ture should  therefore  be  interpreted  in  the  same  sense 
as  choosing  the  atmospheric  ice  pressure  in  the  formula- 
tion of  the  in-situ  freezing  temperature. 


ANALYSIS 

We  shall  use  (4)  of  segregation  freezing  temperature 
as  one  of  the  boundary  conditions  of  the  simultaneous 
flows  of  heat  and  water  to  analyze  the  formation  of  a 
single  ice  lens.  We  shall  make  the  physical  system  as 
simple  as  possible  to  keep  the  analysis  feasible. 

We  assume  the  unfrozen  soil  underlying  the  ice  lens 
to  be  incompressible  under  the  action  of  the  surcharge 
and,  moreover,  under  the  action  of  the  flows  of  heat 
and  water.  At  present  this  assumption  is  needed  be- 
cause the  currently  available  water  flow  equations  do 
not  include  volume  change  caused  by  absorption  and 
depletion  of  water.  Also,  we  do  not  yet  know  the  con- 
stitutive equations  of  soils  to  describe  the  deformation 
due  to  surcharge  and  water  content  variations.  Unifi- 
cation of  hydraulics  and  mechanics  still  seems  to  be  a 
remote  goal.  The  assumption  of  incompressibility 
obviates  these  difficulties.  Furthermore,  this  assumption 
simplifies  the  analysis,  because  the  segregation  freezing 
front  overlying  an  incompressible  unfrozen  soil  layer 
stays  at  the  initial  level  until  in-situ  freezing  replaces  the 
ongoing  segregation  freezing. 

In  this  system  the  freezing  front  starts  to  descend 
when  in-situ  freezing  begins.  The  selection  rule,  suting 
which  of  the  two  processes  should  start,  emerges  at  the 
end  of  the  analysis. 

We  assume  that  segregation  freezing  takes  place  at  the 
ground  surface.  Then,  we  may  not  consider  the  compli- 
cated flow  of  unfrozen  water  in  the  frozen  region.  In 
fact,  we  are  going  to  analyze  the  frost  needle  formation 
on  the  ground  surface. 
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According  to  the  present  theory,  frost  needles  grow 
on  soil  particles.  Pore  water  between  soil  particles  may 
or  may  not  freeze,  because  7,  £ 7|.  In  this  analysis, 
we  disregard  the  individuality  of  the  frost  needles  (as 
seen  on  the  cover)  and  suppose  that  the  ice  lens  forma- 
tion and  heat  and  water  flows  are  uniform  in  the  hori- 
zontal direction.  In  other  words,  we  suppose  that  the 
flows  are  one-dimensional  in  the  vertical  direction.  The 
aim  of  this  analysis  is  not  the  formulation  of  actuality 
but  the  clarifying  of  the  implication  of  our  assumptions. 

We  will  analyze  only  for  the  limit  of  / -►  0.  For 
r 0,  we  can  linearize  the  highly  nonlinear  equations 
of  simultaneous  flows  of  heat  and  water,  and  can  solve 
them  analytically. 

Before  entering  into  the  details  of  the  analysis,  it  is 
appropriate  to  give  an  overview  of  the  analysis. 

First,  we  shall  solve  the  heat  conduction  in  the 
nascent  ice  layer  by  Portnov’s  (1962)  method,  of  which 
the  essence  is  given  in  Appendix  A.  The  boundary 
temperature  conditions  are  the  step-change  air  tem- 
perature Tf^  at  the  upper  side  of  the  ice  \ensAA  in 
Figure  3,  where  X = -/»(r)  and  the  segregation  freezing 
temperature  Tj  at  its  lower  side  55,  where  X = 0.  This 
solution  enables  us  to  express  the  temperature  gradient 
at  55  as  a function  of  7^  and  7,. 

Second,  we  shall  solve  the  unsaturated  water  flow  in 
the  unfrozen  region;  i.e.  we  shall  determine  the  water 
content  W(x,f).  The  boundary  condition  at  x = 0 is 
that  the  water  content  at  X = 0 suddenly  drops  to  a 
certain  unknown  value  W(0,0)  at  the  outset  of  the  ice 
lens  formation.  We  assign  an  arbitrary  number  W(0,0) 
to  the  boundary  value  at  x = 0 and  r = 0.  The  initial 
condition  is  that  W(x,0)  = constant  for  0 < x < <*>. 
The  boundary  condition  at x = “is  that  W(«>,r)  = 
constant.  These  two  constants  must  obviously  be 
equal  to  each  other.  The  solution  of  this  problem  en- 
ables us  to  calculate  the  flux  of  water  entering  the 
freezing  front.  All  this  water  becomes  ice  to  form  frost 
needles;  thus,  we  can  calculate  the  ice  lens  growth  rate 
dhidt. 

Third,  we  shall  solve  the  equation  of  the  double 
heat  transfer,  convected  by  the  water  flow  and  con- 
ducted through  the  soil  mass,  by  using  the  segregation 
freezing  temperature  given  by  (4)  as  one  of  the 
boundary  conditions.  We  evaluate  the  thermal  con- 
ductivity and  heat  content  in  the  duplicate  heat  trans- 
fer equation  by  use  of  the  water  content  distribution 
H'(x,f)  found  above.  The  solution  of  the  duplicate 
heat  transfer  enables  us  to  calculate  the  temperature 
gradient  at  the  segregation  freezing  front  as  a function 
of  the  surcharge  w and  the  boundary  water  content 
value  W(0,0). 
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Figure  3.  Analysis  of  Ice  lens  formation. 


Finally,  we  shall  use  the  energy  balance  equations 
at  the  segregation  freezing  front.  We  substitute  the  ice 
lens  growth  rate  dhidt  and  the  two  temperature  gradient 
equations  previously  formulated  at  both  sides  of  the 
freezing  front  into  the  energy  balance  equation.  Then, 
we  can  find  surcharge  w,  i.e.  the  frost-heaving  pressure, 
in  terms  of  the  air  temperature  7^  and  the  boundary 
water  content  value  W(0,0). 

The  selection  rule  is  given  by  the  expression  of  w. 

If  the  frost-heaving  pressure  w is  zero  or  positive, 
segregation-freezing  begins.  If  the  frost-heaving  pres- 
sure is  negative,  in-situ  freezing  begins.  Soil  data  in 
this  calculation  were  collected  from  many  sources,  and 
the  soils  were  not  incompressible;  however,  the  result 
is  deemed  reasonable. 


Heat  conduction  in  the  ruscent  ice  layer 

The  ice  lens  is  lifted  as  a whole  at  the  rate  of  dhidt. 
Applying  the  theory  of  heat  conduction  in  a moving 
medium  (Carslaw  and  Jaeger  19S9),  we  have  the  dif- 
ferential equation  of  the  temperature  7^  of  the  growing 
ice  lens; 


dh 


327, 


dr  ~ dt  ~5Jr  3x2 
where  Oj  is  the  thermal  diffusivity  of  ice.  Let 
z = X + /7(r) 
and  then  (5)  becomes 


(5) 


37, 

IT 


= Oi 


3z2  ’ 


where  0 S z 5 h{t).  The  solution  of  this  problem  by 
using  Portnov’s  (1962)  method  is  shown  in  Appendix 
A. 
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To  consider  the  limit  of  f -*■  0,  we  take  only  the  first 
(i.e.  n = 0)  terms  in  (A6)  and  (A8)  in  Appendix  A,  and 
approximate  them  by 


'‘'■''="*'■'"‘[1^1 


h(t)  = 2pN/a7. 

where  a,  b,  and  p are  unknown  constants.  To  deter- 
mine a and  b,  we  use  the  conditions: 

= T,  (8 


T-il-ftW.f)  = Tf^,  (9 

where  is  the  air  temperature  on  the  ice/air  inter- 
face AA.  Considering  that 

/»(0)  = 0,  (1 

we  denote  the  value  of  T,  at  A = 0 by 

r,(0)  = r,o-w/(p.,f.n.  (1 

TIics  we  have 

g = r,(0)  - Tf,  tricn 

erfp 


> = (ui 

erfp 

The  constant  p will  be  given  at  the  end  of  the  next 
section. 

Water  flow  in  the  unfrozen  soil 

We  express  the  one-dimensional  flow  of  unsaturated 
water  with  the  following  equation  (Miller  and  Klute 
1967): 

(141 

where  W(x,f)  is  the  degree  of  saturation  by  the  movable 
pore  water  such  that  W = 1 when  the  pore  is  saturated 
with  the  movable  water  and  W = 0 when  no  movable 


water  exists  in  the  pore,  and  D(W)  is  the  hydraulic 
conductivity  of  unsaturated  water  flow  as  a function 
of  W.  We  believe  the  degree  of  saturation  is  more  con- 
venient for  mathematical  analysis  than  the  water  con- 
tent conventionally  used  in  soil  engineering.  The  cross- 
effect of  the  temperature  gradient  on  the  water  flow 
may  not  be  included  in  (14),  because  our  interest  is  in 
the  flow  of  liquid  water,  not  in  the  flow  of  water  vapor 
(Philip  and  DeVries,  1957).  A question,  “What  will 
happen  if  air  is  unavailable  to  the  frost-heaving  system?” 
is  discussed  in  Appendix  B. 

The  flux  v(x,t)  of  liquid  water  is  given  as 

•'(V)  = (1 -Xs)£>(W)(8W/ax)  (15) 

where  Xs  is  the  volume  occupied  by  both  soii  particles 
and  unmovable  water  in  a unit  volume  of  the  soil  mass. 
The  volume  of  movable  water  in  the  unit  voiume  is 
(1  “ Xs)h^-  We  have  formulated  y{x,t)  in  (15)  to  be 
positive  in  the  upward  direction. 

We  assume  the  diffusion  coefficient  D(W)  in  (14), 
following  Gardner  (1958a,  1959),  to  be 

Z)(W)  = Do  exp(^W)  (16) 

where  Dq  and  |3  are  constants.  However,  this  form  of 
D(W)  is  inadequate,  because  it  does  not  become  constant 
in  the  neighborhood  of  W = 1 as  required  by  Darcy  flow 
(i.e.  the  water  flow  saturating  pores).  A few  years  after 
the  analysis  presented  here  was  finished,  we  found 
another  formula  (Gardner  1958b)  that  satisfies  this  re- 
quirement. This  analysis  was  not  revised,  however,  be- 
cause we  believe  that  this  defect  in  our  data  of  water  con- 
tent is  tolerable,  as  will  be  shown  in  the  numerical 
analysis  section.  A more  realistic  analysis  should  be 
attempted  when  the  rigid  soil  frost-heaving  test  will  be- 
come available. 

Use  of  u defined  by 

u = expipW)  (17) 

simplifies  eq  14,  where  D(W)  is  given  by  (16),  to 

^ bu  - b^u  /lor 

In  view  of  Portnov’s  (1962)  formulation  (App.  A), 
we  may  assume  u{x,t)  to  be  in  the  following  form: 


u{x.t)  = x;  (DoD'-ZZ  D„(t) 


where  f is  defined  by 


i = x/y/2D^  (20) 

and  U„{i)  is  a function  of  f only.  To  discuss  the  limit 
t -*  0,  we  need  only  the  first  term: 


uM  = Uoi^). 


(21) 


Then  (18)  becomes 


dUo 

di 


d^Ufi 

dk^ 


= 0. 


(22) 


Use  of  a single  independent  variable  { demands  that  the 
boundary  condition  at x = “at  time  t > 0 and  the 
initial  condition  in  the  region  0 «:  x < “>  reduce  to 
a single  condition: 


^oM  = exp(/3W'„) 


(23) 


of  Repeating  the  analytical  continuation  as  many 
times  as  necessary,  one  can  find  the  value  Uq(“), 
which,  however,  in  general  is  different  from  the  given 
value  of  Uq(°°).  One  renews  t/Q(0)  and  repeats  the 
same  procedure  until  Uq(‘»)  agrees  with  the  given  value. 
Thus,  the  Scott  and  Hanks  method  enables  us  to  express 
t/o(0)  numerically  in  terms  of  Uo(0).  In  this  way  we 
can  find  the  numerical  solution  of  Uq{^)  involving 
1V(0,0)  as  a parameter. 

We  may  approximate  ^(x,^)  by  substituting  in 
(25)  for  Win  (15); thus,  we  find 

vix  t)  = (1  - X ) -J—  (29) 

by  use  of  (26).  The  balance  of  mass  at  the  freezing 
front  X = 0 is  given  by 

= {3«') 


where  W„  is  a constant  such  that 

•♦^(0.0)  < W_  < 1.  (24) 

We  determine  the  boundary  value  Uq(0)  as  follows. 
Use  of  (19)  in  (1 7)  shows  us  that  W(x,f)  can  be  ex- 
pressed with  a series  similar  to  (19).  For  the  limit  of 
/ -►  0,  taking  only  the  first  term,  W(x,f)  may  be 


approximated  by 

=.  (25) 

where  v^od)  is  related  to  by 

i^o(«)  = exp(0.^o(«))-  (26) 

Because  ^o(^)  continuous,  it  must  satisfy 

W(0,0)  = ^|>o{0).  (27) 

Thus  we  get 

Uo(0)  = exp(^W(0,0)).  (28) 


To  solve  (22)  with  the  boundary  conditions  (23)  and 
(28),  we  used  Scott  and  Hank’s  (1962)  method.  Assuming 
an  arbitrary  value  of  Uo(0)  and  given  Uo(0)and 
£7j)(0),  one  can  compute  the  higher  derivatives  U^HO) 

(n  > 2)  by  use  of  (22).  Thus,  one  can  formulate  a 
Taylor  series  in  the  neighborhood  of  { =0.  If  the 
convergence  deteriorates  as  | increases,  one  can 
analytically  continue  the  Taylor  series  to  a new  series 
that  better  converges  in  a new  range  of  larger  values 


Use  of  (29)  in  (30)  yields  the  differential  equation  of 
h(t),  which  on  integration  with  the  initial  condition 
/;(0)  = 0 yields  p introduced  in  (7): 

= (1  - Xs)  ^ <^(0).  (31 ) 

Heat  transfer  in  the  unfrozen  soil 

We  shall  formulate  the  equation  of  the  double 
heat  transfer,  converted  by  the  water  flow  and  con- 
ducted through  the  soil  mass.  The  heat  content  Cj  per 
unit  volume  of  unfrozen  soil  mass  is 

C2  = C^W(1-Xs)+c,x,  (32) 

where  Cj  is  the  heat  content  per  unit  volume  of  soil 
particles  including  unmovable  water,  and  the  heat 
content  per  unit  volume  of  water.  Let  T2(x,r)  be  the 
temperature  of  the  unfrozen  soil  mass.  The  convective 
heat  flow  formulated  by 

Q = (33) 

is  positive  upward  because  the  flux  v(x,r),  defined  by 

(1 5),  is  positive  upward.  The  conductive  heat  flow  R 
is  given  by 

R = (34) 

where  K2  is  the  thermal  conductivity  of  the  soil  mass. 
Because  X is  positive  downward,  ^ is  positive  down- 
ward. The  duplicate  heat  transfer  is  formulated  by 
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31c2  T2)/3t  = aO/3x  - bRIbx. 

To  simplify  (35),  we  derive  the  relation 


(35) 


dc2ldt  = dv/dx  (36) 

by  differentiating  (32)  with  regard  to  t and  using  (14) 
and  (1 5)  on  the  assumption  thatc,  x$  = constant,  i.e. 
that  soil  particles  including  unmovable  water  do  not 
move.  Substituting  (33)  and  (34)  into  (35)  and  using 
(36),  we  find 


dT2 

W 


dT2 

’lx" 


(37) 


It  may  be  verified  that  the  equation  of  the  triple 
heat  transfer  (Philip  and  DeVries  1957,  and  DeVries 
1958),  containing  convection  by  water  vapor,  con- 
vection by  liquid  water,  and  conduction  through  the 
soil  mass,  reduces  to  (37)  when  the  former  is  simplified 
on  the  assumption  that  the  flow  of  water  vapor  is 
negligible. 

Equation  (37)  of  the  duplicate  heat  transfer  includes 
the  effect  of  the  water  flow  v on  the  temperature 
gradient,  but  (14)  of  the  unsaturated  water  flow  does 
not  include  the  reciprocal  relationship,  i.e.  the  effect 
of  the  temperature  gradient  on  the  water  flow.  The 
validity  of  Onsager’s  reciprocal  relationship  (DeGroot 
and  Mazur  1962)  is  not  claimed  in  this  paper;  nor  does 
this  relationship  hold  in  the  theory  by  Philip  and 
DeVries  (1957)  and  DeVries  (1958).  The  neglect  of 
the  relationship  seems  to  be  natural,  although  not  yet 
proven,  in  the  simultaneous  flows  of  heat  and  water 
through  soil,  because  heat  can  penetrate  soil  particles 
but  water  cannot  — a condition  that  is  not  considered 
in  theoretical  physics  for  proving  Onsager’s  relation- 
ship. 

We  cast  (32)  into  a form  convenient  to  the  heat 
flow  analysis 

C2  = (1  -r*rW)  (38) 

by  introducing  two  constants  q and  r through  the  fol- 
lowing two  equations. 


qr  = ^ - Xs 

q(l  -r)  = c^xjc^. 


(39) 


whose  right-hand  side  is  the  one  found  by  letting 
W = 1 in  the  right-hand  side  of  (32).  The  constant  r 
is  in  the  range  0 < r < 1 , because  dividing  the  first 
equation  of  (39)  with  the  one  at  the  top  of  this  column 
shows  that/-  = c^(l  -Xs)/(c^(1  - Xs)  "'■c,  Xsl  • 

We  used  Kersten’s  (1949)  equation 

K2  = *20  ^0  '08 {^) 

to  express  the  thermal  conductivity  K2  of  unsaturated 
soil.  The  constant  X20  is  the  thermal  conductivity  for 
the  saturated  condition  (W  = 1 ),  given  by 

*20  = Pw  <7®20  (41) 

where  020  is  the  thermal  diffusivity  of  the  saturated 
soil.  The  constant  Xg  is  a soil  constant. 

In  view  of  Portnov’s  formulation  (App.  A),  we  may 
assume  T{x,t)  to  be  in  the  following  form; 

OP 

T2M  = (42) 

n=0 

where  we  have  introduced  a function  of  | only. 
Taking  the  lowest  term  of  t,  we  approximate  T2{x,t) 
with 

T2M  k eo(f).  (43) 


Substituting  T2  from  (43),  vfrom  (29),  C2  from  (38), 
and  K2  from  (40),  (37)  transforms  to 

<»20  tf^Q0^L^).°2O  L 

^0  L ^ ^0  1 A IOgWo({)J 


X ^ = 0 (44) 


where  we  have  defined 


m 


1-^XIogWo(£) 


(45) 


The  boundary  condition  at  x = 0 is 


00  (0)  = r,(0). 


r 

I 

i' 


The  constant  qc^  expresses  the  heat  contained  in  the 
water-saturated  soil  mass.  To  prove  this,  note  that  the 
sum  of  the  two  equations  in  (39)  yields  the  relation 
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(46) 


(54) 


The  boundary  condition  at  x = <»  and  the  initial  con-  T, 

dition  for  0 < X < » reduces  to  a single  condition  p.  L 


= r.(o). 

Equation  44  integrates  to 

t 

= ©o(0)(1  +X  logM'o(O))  / 
0 


G(t?)  dv  + 

+ UO) 


(47) 


(48) 


/ G(7j)r/Tj 

_a 


l! e ' 


-u2 


. erf#r 


V25^  / G(-n)<i-n 


where  0o(O)  is  the  yet  unknown  value  of  dQQ{()/di 
at  { = 0,  and  G{ri)  is  defined  by 

The  boundary  condition  (46)  is  satisfied  by  (48).  The 
boundary  condition  (47)  is  satisfied  if  9j)(0)  is  chosen 
to  be 


Energy  balance  at  the  segregation-freezing  front 
The  balance  of  energy  at  the  segregation  freezing 
front  is  described  by 


where  (*2)0  *2  tout'd  from  (40)  by  let- 

ting = M'(0,0).  Notations  (3r,/3x)Q  and  (97'2/3x)q 
are  values  of  3r,/dx  and  BTj/dx  at  X = 0.  For  the 
limit  of  Z -♦  0 we  find 


l^Jl)  - ^(0)-^  1 

' 3^  / erfp 


(52) 


by  substituting  (52)  and  (53)  for  the  temperature  gradi- 
ents in  (51 ),  (7)  for  />(f),  and  (1 1 ) for  T^(0)  in  (52).  In 
this  equation  p is  a function  of  11^(0, 0)  that  can  be  found 
by  expressing  t/J)(0)  in  (31 ) as  a function  of  1^(0, 0), 
which  we  have  derived  by  use  of  Scott  and  Hanks’ 
method;  Tj  is  also  a function  of  1^(0, 0)  as  given  in  (55) 
below. 

Equation  (54)  gives  the  selection  rule.  If  the  right- 
hand  side  of  (54)  is  zero  or  positive,  segregation  freezing 
starts.  If  it  is  negative,  segregation  freezing  cannot  start 
but,  instead,  in-situ  freezing  begins. 

Numerical  computation 

The  only  experiment  on  segregation  freezing  in- 
cluding the  measurement  of  unsaturated  water  flow  is, 
to  our  knowledge,  Hoekstra’s  (1966, 1967).  The  soil 
he  used  was  Fairbanks  silt,  which  does  not  satisfy  the 
assumption  of  rigid  pores;  therefore  (4)  of  may  not 
be  exact  for  this  soil.  The  water  flow  was  unsaturated, 
the  initial  W being  equal  to  0.82;  therefore  (16)  may  be 
used  for  D(W).  Although  his  soil  column  was  of  finite 
length,  we  may  use  his  data  in  our  analysis,  because  we 
consider  only  the  limit  oft  -*  0.  The  porosity  was 
0.36  and  therefore  Xj  = 0.64.  We  determined  the 
specific  density  Pj  of  the  soil  by  equating  Pj  Xs  to  the 
dry  density,  which  was  1670  kg/m^. 

Low  et  al.  (1968)  observed  almost  complete  linearity 
between  Tf  and  W for  a Wyoming  Na-bentonite  nearly 
saturated  witii  water.  Assuming  that  this  relationship 
is  valid  even  for  other  soils,  we  formulated 


from  (6)  by  use  of  (13),  and 


r,  = To-vd-W) 


(55) 


from  (43)  by  use  of  (50).  Thus,  we  can  express  w in 
terms  of  parameters  W(0,0),  T.,  and  as 


where  Tq  = 273.1 5 K and  v is  a soil  constant.  We 
chose  V = 0.140,  referring  to  Keune  and  Hoekstra 
(1967). 

Kersten’s  (1949)  data  of  the  thermal  conductivity  of 
Fairbanks  silt  gave  = 0.892  and  K20  = 1.6039  W/mK. 
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Figure  4.  Frost-heaving  pressure  w as  function 
of  and  W(0,0). 


Table  I.  Temperatures  and  degree  of  saturation  for 
needle  ice  formation. 


W(0,0) 

T^rc) 

r,  " r/o 

0.4 

-0.193 

-0.084 

0.5 

-0.156 

-0.07 

0.6 

-0.115 

-0.056 

0.7 

-0.0732 

-0.042 

0.8 

-0.0347 

-0.028 

0.82 

-0.0279 

-0.025 

The  specific  heat  of  the  dry  soil  given  by  him  was 
79S  J/kg  K;  therefore,  c,  is  given  by  c,  = 795  p,. 

Using  (39)  we  find  q = 0.6853  and  r = 0.5370. 

Values  of  Dq  and  0 determined  by  Hoekstra’s  (1%6, 
1967)  data  were  Og  = 2.92x10"*  m^/s  and  0 = 2.88. 
The  method  of  determination  of  these  values  is  not 
mentioned  here,  because  special  knowledge  of  un- 
saturated water  flow  was  used  for  their  determination 
as  mentioned  by  us  (1970).  The  values  of  Dg  and  0 
were  reasonable  as  compared  with  other  soils. 

Numerical  computation  was  performed  keeping 
T".  = 5‘’C  constant  and  varying  and  M'(0,0).  The 
relation  between  w and  is  shown  in  Figure  4 with 


ll'(0,0)  as  a parameter.  The  ice  lens  that  forms  when 
w = 0,  i.e.  under  atmospheric  pressure,  is  usually 
called  needle  ice.  Under  this  condition,  T,  is  equal  to 
r,,  which  may  be  computed  from  (55)  by  using  M'(0,0) 
for  IF.  The  values  thus  found  are  shown  in  Table  I and 
Figure  4. 

Although  we  have  used  the  simplifying  assumption 
of  rigid  pores  and  collected  the  input  data  from  a 
variety  of  sources,  the  results  shown  in  Figure  4 and 
Table  I are  reasonable  when  compared  with  observa- 
tions in  the  laboratory  and  in  nature. 


LITERATURE  CITED 

Bell,  E.T.  (1934)  Exponential  polynomials.  Annals  of 
Mathematics,  vol.  35,  no.  2,  p.  258-277. 

Carslaw,  H.S.  and  ).C.  Jaeger  (1959)  Conduction  of  heat  In 
solids.  2nd  ed.,  Oxford:  Clarendon  Press. 

Corte,  A.E.  (1962)  Vertical  migration  of  particles  in  front  of 
a moving  freezing  plane.  Journal  of  Geological  Research, 
vol.  67,  no.  3,  p.  1085-1090. 

OeGroot,  S.R.  and  P.  Mazur  (1962)  Non-equillbrium  thermo- 
dynamics. New  York:  Interscience  Publishers. 

DeVries,  O.A.  (1958)  Simultaneous  transfer  of  heat  and  mois- 
ture in  porous  media.  Transactions  of  American  Geo- 
physical Union,  vol.  39,  no.  5,  p.  909-916. 

Everett,  D.H.  and  ).M.  Haynes  (1965)  Capillary  properties  of 
some  model  pore  systems  with  special  reference  to  frost 
damage.  RILEM  Bulletin  (Reunion  Internationale  des 
laboratoires  d'euais  et  de  recherches  sur  let  matiriaux  et 
les  constructions).  New  Series,  no.  27,  p.  31-38. 

Gardner,  W.R.  (1958a)  Mathematics  of  isothermal  water  con- 
duction in  unsaturated  soil.  Highway  Research  Board 
Special  Report  40,  p.  79-87. 

Gardner,  W.R.  (1958b)  Some  steady-state  solutions  of  the 
unsaturated  moisture  flow  equation  with  application  to 
evaporation  from  a water  table.  Soil  Science,  vol.  85, 
p.  228-232. 

Gardner,  W.R.  (1959)  Solutions  of  the  flow  equation  for  the 
drying  of  soils  and  other  porous  media.  Soil  Science  Society 
of  America  Proceedings,  vol.  23,  p.  183-187. 

Gautschi,  W.  (1964)  Error  functions  and  Fresnel  integrals.  In 
Handbook  of  Mathematical  Functions,  (M.  Abramowitz 
and  I.  Stegun,  Eds.)  National  Bureau  of  Standards  Applied 
Mathematics  Series,  55.  Washington,  O.C.,  United  Sutes 
Department  of  Commerce,  p.  295-329. 

Glen,  ).W.  (1974)  The  physics  of  ice.  CRREL  Cold  Regions 
Science  and  Engineering  Monograph  ll■C2a.  AO  778009. 

Hobbs,  P.V.  (1974) /cepA>'j<cs.  London;  Oxford  University 
Press. 

Hoekstra,  P.  (1966)  Moisture  movement  in  soils  under  tempera- 
ture gradients  with  the  cold  side  temperature  below 
freezing.  Water  Resources  Research,  vol.  2,  no.  2,  p.  241- 
250. 

Hoekstra,  P.  (1967)  Moisture  mover,-  ent  to  a freezing  front. 

In  Geochemistry,  Precipitation  *.vaporatlon,  Soll- 
Molsture,  Hydrometry,  Interna.,  nal  Association  of 
Scientific  Hydrology  (Bern),  pub.  78,  p.  411-417. 

Jackson,  F.  (1964)  The  solution  of  problems  involving  the 
melting  arid  freezing  of  finite  slabs  by  a method  due  to 


9 


i 

I 


\ 


ii 


[ 

[ 

[ 

\ 

! 

F 

r 


Portnov.  Proceedings  of  the  Edinburg  Mathematical 
Society,  vol.  14,  set.  II,  part  2,  p.  109-128. 

Kersten,  M.S.  (1949)  Laboratory  research  for  determination 
of  the  thermal  properties  of  soils.  Final  Report  to  the 
Corps  of  Engineers,  U.S.  Army  District,  St.  Paul,  Min- 
nesota. 

Keune,  R.  and  P.  Hoekstra  (1967)  Calculating  the  amount  of 
unfrozen  water  in  frozen  ground  from  moisture  character- 
istic curves.  CRREL  Special  Report  1 14.  AO  659781. 

Koopmans,  R.W.R.  and  R.O.  Miller  (1966)  Soil  freezing  and 
soil  water  characteristic  curves.  Soil  Science  Society  of 
America  Proceedings,  vol.  30,  p.  680-685. 

Loch,  J.P.G.  and  R.O.  Miller  (1975)  Tests  of  the  concept  of 
secondary  frost  heaving.  Soil  Science  Society  of  America 
Proceedings,  vol.  39,  p.  1036-1041. 

Low,  P.F.,  D.M.  Anderson  and  P.  Hoekstra  (1968)  Some 
thermodynamic  relationships  for  soils  at  or  below  the 
freezing  point.  1:  Freezing  point  depression  and  heat 
capacity.  Water  Resources  Research,  vol.  4,  no.  2, 
p.  379-394. 

Miller,  E.E.  and  A.  Klule,  (1967)  The  dynamics  of  soil  water. 

1:  Mechanical  forces.  (\n  irrigation  of  Agricultural  Lands, 
ch.  13.)  American  Society  of  Agronomy,  Madison, 
Wisconsin. 

Penner,  E.  (1967)  Heaving  pressure  in  soils  during  unidirectional 
freezing.  Canadian  Geotechnical  journal,  vol.  4,  p.  398-408. 

Philip,  ).P.  and  D.A.  DeVries  (1957)  Moisture  movement  in 
porous  materials  under  temperature  gradients.  Transactions 
of  American  Geophysical  Union,  vol.  38,  p.  222-232. 

Portnov,  I.G.  (1962)  Exact  solution  of  freezing  problem  with 
arbitrary  temperature  variation  on  fixed  boundary. 

Soviet  Physics-Doklady,  vol.  7,  no.  3,  p.  186-189. 

Riordan,  |.  (1946)  Derivative  of  composite  functions.  Bulletin 
of  the  American  Mathematical  Society,  vol.  52,  p.  664-667. 

Riordan,  J.  (1949)  Inversion  formulas  in  normal  variable 
mapping.  Annalen  of  Mathematical  Statics,  vol.  29,  p.  419. 

Schofield,  R.K.  (1935)  The  pF  of  the  water  In  soil.  Trans- 
actions of  the  3rd  International  Congress  of  Soil  Science, 
vol.  2,  p.  37-48. 

Scott,  E.).  and  R.).  Hanks  (1962)  Solution  of  the  one- 
dimensional diffusion  equation  for  exponential  and  linear 
diffusivity  functions  by  power  series  applied  to  moisture 
flow  in  soils.  Soil  Science,  vol.  94,  no.  5,  p.  314-322. 

Sutherland,  H.B.  and  P.N.  Gaskin  (1973)  Pore  water  and 
heaving  pressures  developed  in  partially  frozen  soils. 
Permafrost,  2nd  International  Conference,  North 
American  Contribution,  National  Academy  of  Sciences, 
Washington,  D.C.,  p.  409-419. 

Takagi,  S.  (1959)  Theory  of  freezing-point  depression  of  soil 
water,  and  a note  on  the  extra-thermodynamics  of  soil 
moisture.  Soil  Science,  vol.  88,  no.  1,  p.  25-31. 

Takagi,  S.  (1963)  Fundamentals  of  the  theory  of  frost-heaving. 
Permafrost,  1st  International  Conference,  National  Academy 
of  Sciences-National  Research  Council,  Washington,  D.C., 
p.  203. 

Takagi,  S.  (1965)  Principles  of  frost  heaving.  CRREL  Re- 
search Report  140.  AD  626174. 

Takagi,  S.  (1970)  An  analysis  of  Ice  lens  formation.  Water 
Resources  Research,  vol.  6,  no.  3,  p.  736-749. 

Takagi,  S.  (1974)  Heaving  of  Fairbanks  slit  during  freezing 
with  and  without  air  available.  CRREL  Technical  Note 
(unpublithed). 

Takagi,  S.  (1975)  Fundamentals  of  ke  lens  formation.  15th 
National  Heat  Transfer  Conference,  San  Francisco, 
California,  Paper  no.  28. 


Takagi,  S.  (1977)  Segregation-freezing  temperature  as  the 
cause  of  suction  force.  Proceedings  of  the  International 
Symposium  on  Frost  Action  in  Soils,  Lulea,  Sweden, 
vol.  1,  p.  59-66. 

Terzaghi,  K.  (1942)  Theoretical  soU  mechanics.  New  York: 
)ohn  Wiley  and  Sons. 

Williams,  P.).  (1964)  Unfrozen  water  content  of  frozen  soils 
and  toils  moisture  suction.  Geotechnique,  vol.  14,  no.  3, 
p.  231-246. 


10 


APPENDIX  A 


f 


i 

r 

k 


Essence  of  Portnov’s  method 

Portnov  (1962)  presented  an  interesting  idea  for 
handling  the  moving  boundary  heat  conduction 
problems.  Jackson  (1964)  extended  it  and  showed 
several  examples.  However,  their  presentations  are 
still  complicated.  In  the  following  a simple  formula- 
tion of  Portnov’s  idea  is  presented. 

It  is  well-known  that  a solution  of  the  heat  con- 
duction equation 


il  - 

dt 


a 


3x2 


in  the  infinite  region,  - “>  < x < oo,  is 


' vir  / I''" 

where  0(x)  represents  the  initial  temperature  distribu- 
tion. Use  of  (A1 ) enables  one  to  find  the  solution  in 
the  moving  boundary  region  0 ^ x S h{t). 

Let  0(x)  = 0]  (x)  in  0 < X < > and  0(x)  = 
02{x)in-«>  < X < 0.  Assuming  that  0,  (x)  and 
02  (x)  are  analytic  in  their  respective  regions,  let 


where  /('’)erfcx  is  the  n-time  repeated  integral  of  erfcx, 
which  may  be  expressed  in  the  form  of  a single  integral; 

/<")erfcx  = f e-"^  du.  (A3) 

^ J n\ 

To  find  (A2),  let  f - x = 2u>/od  in  the  region 
0 < f < o®  and  X - f = 2u\/ad  in  the  region 
- “>  < t < 0.  Then  use  of  (A3)  easily  yields  (A2). 

The  customary  notation  /""erfcx  is  rejected  here  in 
favor  of  /^'’ferfcx,  because  in  the  functional  nota- 
tion can  be  confused  with  (\/^)‘". 

The  formula  changing  the  negative  argument  of 
/<">erfc(-x)  to  the  positive  argument  of  /("lerfcx 
(c.f.  formula  7.2.11  of  Gautschi  1964)  simplifies  to 

(-  1)"  /("(erfcx  + /('’lerfc(-x)  = — 3 — EM),  (A4) 

2"-’/>! 

when  the  polynomial  ^^(x),  defined  by 

£„{x)  = e-^^  (AS) 


01 W =2]  0i”>x" 

n=0 

•• 

<>2^)  = 2^  0^>x". 
n=0 

Then,  one  can  integrate  (Al)  to 

•• 

T{^,t)  = j 2^  /»!  {2>/caY'^ip\'’)  /("lerfc^- + 

(A2) 


= l-'’H„{/x) 

is  introduced,  where  H„{x)  is  the  Hermite  polynomial 
and  / = \/^.  Using  (A4),  one  can  transform  (A2)  to 
contain  positive  arguments  only: 

n...i  ■ 

Equation  (A6)  shows  that  the  temperature  functions 
for  X = 0 and  x = h{t)  are 
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(A7) 


«• 

mo  = E 

n=0 

nmo  = 

n=0 

Therefore,  (A6)  can  express  the  temperature  in  the 
growing  ice  lens,  if  />(0  is  a power  series  of  y/F : 

m 

h(0  = , (A8) 

n=l 

The  series  must  begin  with  y/T  because  is  in  the 
arguments  of  functions  5„(x/(2v^ ))  and 
/('')crfc(x/{2\/oF)).  Bell’s  formula  (Bell  1934,  Riordan 
1946, 1949)  which  gives  the  nth  derivative  of  a func- 
tion %(x)),  may  be  used  to  express  a„  and  b„  in 
(A7)  in  terms  of  the  derivatives  of  E^(h(t)l(2\/at)) 
and  /('’)erfc(/»(r)/(2N/^tf ))  and  arbitrary  constants 
0^”)  and 
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APPENDIX  B 

Frost-heaving  without  air  available 

We  can  theoretically  prove  that  segregation  freezing 
cannot  start  in  a rigid  soil  whose  pores  are  saturated 
with  de-aired  water. 

To  prove  this,  we  may  assume  that  the  flow  is  one- 
dimensional. The  equation  of  continuity,  div  V = 0,  re- 
duces in  one-dimensional  flow  to  bvjbx  = 0,  where  V is 
the  velocity  vector,  v the  vertical  component,  and  x the  I 

vertical  coordinate.  Therefore,  v is  a function  of  t only. 

Portnov’s  formulation  shows  that  ice  thickness  h for 
an  initial  small  period  is  proportional  to  y/F.  Sub- 
stitution of  this  result  into  (30)  of  the  balance  of 
mass  at  the  freezing  front  indicates  that  v for  the 
initial  small  period  must  be  proportional  to  r^.  In 
other  words,  initial  velocity  is  infinite  throughout 
the  entire  domain,  0 k x > <».  Under  this  initial 
condition,  the  problem  of  water  flow  cannot  be  solved 
and,  therefore,  segregation  freezing  cannot  start. 

We  showed  experimenully  (Takagi  1 974)  that  this 
theoretical  conclusion  with  regard  to  rigid  soils  does 
not  necessarily  hold  true  with  regard  to  deformable 
soils.  However,  the  mechanics  of  water  flow  in  de- 
formable soils  is  not  readily  understandable. 
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